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Abstract 

We analyze the effects of the space-like penguin diagrams in two body B decays 
B~ — > K~K° and B~ — > tt~K°. Special attention is paid to the operator Qq, whose 
contribution is large in the BSW model, but vanishes in the PQCD method in the ap- 
proximation of neglecting the masses in the final states. Suppressions of the space-like 
penguin diagrams are found in the PQCD method, which implies that the contributions 
of the space-like penguins are small compared to those of the time-like penguins. 
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Penguin diagrams in exclusive B decays have been of great interests in the study of CP 
violations 0. Up to now, most of the quantitative investigations concentrate on the so-called 
time-like penguin diagrams, except in [0, |3| where, using the BSW modeljD, the space-like 
penguin diagrams are found to be important. 

An important feature of the BSW model in nonleptonic decays is the use of the factor- 
ization hypothesis [Q. Under this hypothesis, one decomposes the decay amplitude into the 
product of two hadronic matrix elements which are related to the leptonic or semileptonic 
processes. While this hypothesis works for the tree level processes if some phenomenological 
parameters are introduced to fit the experimental data^, |5|, there exists no firm foundation 
to use it in processes involving penguin diagrams. Because the penguin operators have dif- 
ferent Lorentz structures from the tree level (V — A) x (V — A) ones, it is necessary to use 
the equations of motion for the free quarks to relate the hadronic matrix elements of these 
two kinds of operators^. However, this procedure is not proved since these equations must 
be used for the quarks inside hadrons. Consequently, for a special operator Qe[F]]> the decay 
amplitudes of the space- like penguins always involve the factor M^/(m 9 — m q t), where q and 
q' are light quarks 0, 0, [|. This is a huge factor and more seriously, it might produce the 
singularity at m q = m q i which cannot be cured by the BSW model itself. 

In this work we use the perturbative QCD(PQCD) method]!, || to re-analyze the pen- 
guin dominated processes B~ — > K~K° and B~ — > tt~K°. The calculations of the decay 
amplitudes between the involved hadrons are replaced by those of the scattering amplitudes 
between the quarks, and the non-perturbative aspects of QCD are attributed to the wave 
functions which are taken as external inputs. Since there exists no reliable way to determine 
these wave functions up to now, some approximate forms (like those to be used below) were 
chosen in the literature as practical ways of performing phenomenological analysis. Hence 
the results are model dependent, like those from the BSW model[f|]. 

Following ||, in the PQCD method there is no need of using the factorization hypothesis 
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or the equation of motion in the nonleptonic processes. We calculate the decay amplitudes 
for the penguin diagrams, and compare the contributions from the space-like and time-like 
penguins. We find that within PQCD there is no singularity which exists in the BSW model, 
and the space-like penguins are generally suppressed compared to the time-like ones. 



We take the pure penguin processes as the examples in the following. To start, we use the 
effective Hamiltonian for the charmless nonleptonic B decays as 

G ^ 

Keff = —TkYI V CKMCi(m b )Qi, (1) 

V 2 i=l 

where Vckm denotes the CKM factors, and the operators are 

Qi = {q a up) v ^ A {upb a ) v _ A , Ql = (qu) V - A (ub) v _ A , 
Qs= (qb) V - A ^2(q'q') V - A , Q 4 = (q a b/ 3 ) v „ A J2% ( l' a )v-A, 



<i <i 
Qs= (qb) V - A J2(Q'<l')v+A, Qe= (q a h)v-A^2%q' a )v+A 

q> q> 



(2) 



for q = d, s and q' = u, d, s. The operators Q% to Qq come from the penguin diagrams. 
Performing the Fierz transformations, we get 

Ql = (uu) V - A (qb) v „ A , Ql = (qu) v _ A (ub) v ^ A , 

Qs = J2(Q b )v-A(q'q')v-A, Qi = J2(y' b )v-A(qq')v-A, (3) 

q> q> K > 

Qb= ^2(qb) V -A(q'q')v+A, Q&= -2j2(v' b )s-p(qq')s+p- 

q' q> 

The Wilson coefficients of ([[]) are|7j 

d(m b ) = -0.184, C 2 (m b ) = 1.078, 

C 3 (m b ) = 0.013, C 4 (m 6 ) = -0.035, (4) 
C 5 (m b ) = 0.009, C e (m b ) = -0.041. 
Then the amplitudes for the 2-body nonleptonic B decays are 

< YZ\H eff \B >= ^ E V CKM Ci(m b ) < YZ\Q t \B > . (5) 

v 1 i=l 



In processes with penguin contributions, the amplitudes of @ receive contributions from 
both time-like and space-like penguin diagrams, which are depicted in Figure 1 (a) and 1 
(b). 



Working in the BSW model, the treatments of the space-like penguin contributions of the 
operator Q e need to be discussed in some details. These treatments can be described as the 
following. We use the basis in (^). First, the factorization hypothesis is made to get 

< Y(q'q v )Z(q v q')\Q 6 \B(bq') >= -2 < Y(q'q v )Z(q v q')\(qq') s+P \0 X 0| (q'b) s _ P \B(bq') > . 

(6) 

Second, use of the equations of motion || [| || 

_ , -id^(q^q') _, -id^(q%j 5 b) 

QQ = — , Ql5b= — - (7) 

is needed to relate the two hadronic matrix elements on the right-handed- side of (§) to those 
with the usual ones with V±A currents. Third, the quark currents on the right-handed- sides 
of (0) are approximated by their hadronic counterparts as|4[] 

qb = (qb) H , (8) 

where the sublabels denote the hadronization. And fourth, partial differential operations on 
the hadronic matrix elements are used giving the final result 

2m 2 

<Y(q'q v )Z(q v q')\Q 6 \B(bq')> - 



(m q - m q i)(m b + m q >) 
< Y(q'q v )Z(q v q')\(qq') v \0 >< 0\(q'b) A \B(bq') > . 

(9) 

There exist another ambiguity in the procedure described above besides the usefulness 
of the factorization hypothesis. In using the equations of motion in (0) we have worked at the 
quark level, thus the quark masses appearing in (M) seem to be the current ones. However, it 
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is also allowed to reverse the second and the third steps in the above procedure, which gives 
no explanation of which mass should be used. Furthermore, when q' and q are the same 
flavors then (|9]) results in infinity which violates the analytic property of the amplitude. 
Even without the singularity, the factor m B j{m q — m q i) in (^) enhances the amplitude of Qq 
greatly, which is difficult to understand. 



Now we turn to the PQCD calculations of the penguin amplitudes. In this method, no 
explicit use of the factorization hypothesis is needed. We cannot resolve the ambiguity of 
the quark masses described above. However, this ambiguity is less serious than that in the 
BSW model. 

Following ||, we take the interpolating fields for the mesons as 

*Pb = ^ ^- <j) B (x)~f 5 (tf- m B ), 

*Py = r^^Mi/hs ( 10 ) 

Here I c is an identity in the color space, and we have neglected the masses of the final states. 

The decay amplitudes can be calculated using Figure 2 and 3 for the time-like and the 
space-like penguins, respectively. We neglect terms proportional to the final state masses, 
and keep in each diagram the leading term in the expansion in x±, the momentum fraction 
carried by the light anti-quark in the B meson. Denoting 

< YZ\Qi\B >= A ia + A ib + ■ ■ ■ + A if , (11) 

where A^{k = a,b, ■ ■ ■ f) is the contribution from the k-th diagram of Figure 2 or 3. When 
the diagrams have the topology which can induce the processes, the results are summarized 
in the following, 
(a) The time-like penguins: 



Atr = -|4 J\dx][dy][dz]Mx)Mv)Mz) 4( 2xi + yi + 1) 



3^6 Jo Xi(x! - y 1 )(2x 1 - y ± )' 

A u r + A u r = f\dx][dy][dz]Mx)My)Mz) T 7Z Arr- ~y (12) 

3 V6 Jo xi [xi - yi ) [xi - zi ) 

a time Atime a time r\ 

A 3a — -^36 — ^3/ — u ; 

- 4£ '" = -Si' W[dyl[dzlfe(l) ^ fa)fe(z) ^S)ii^ - 

Atime Atime Atime Atime Atime r\ ("i o\ 



a time a time Atime Atime , 4 time Atime r\ (TA\ 

A 5a — ^56 — A bc — A bd ~+~ A be ~ A 5f ~ U > 

Atime Atime Atime Atime Atime Atime r\ /ir\ 

A 6a — ^66 — A 6c ~ A Qd ~ A 6e ~ A 6f ~ U ? l i0 J 

where [dx] = d£ida;2, [dy] = ck/ick/2, [dz] = dz\dz 2 . 
(b) The space-like penguins: 

a space a space a space . a space ^ 

A 3a — A 3b — A 3e + A Zf ~ U > 

AZ ace = M l\dx}[d y }[dz}Mx)My)<Pz(z)- - -i 

3V6 Jo (1 - yi)(xi - xiyi + xizi - y x z x + y x - 1) 

^ aCe = |4 /Wdy][ds]M^Mv)M*)-fi i 

3V6 Jo zi(l - - zi) 

(16) 



< aCe = f\dx][dy}[dz]Mx)My)Mz)- 7 T^— y 
3V6Jo zi(l-yi) 

= -^/jdx][dy][dz]0 B (x)0 y (y)0 Z (^)-^— y, (17) 

= (k = c,d,e,f), 

Af k ace = (k = a,b, ■■■/), (18) 

^r e = (fc = a,6, ■■■/). (19) 
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In ((L21-J19D g s = ^Aira s is the coupling constant of QCD. 

Independent of the special choices of the wavef unctions, we can draw some general 
conclusions in the PQCD method. First, from (|l9f) , the space-like penguin contribution 
from the operator Q 6 vanishes, and thus the singularity due to the degeneration of the quark 
masses disappears. This is in sharp contrast to the observation within the BSW model 
discussed above. 

Second, comparing the right-hand-sides of equations (|16| - |T9l) and those of equations (|T^ - 
[T5| ) , it is obvious that the contributions from the space-like penguins are always suppressed 
compared to those from the time- like penguins by a factor of x\, the momentum fraction 
carried by the light quark inside the B meson, which is taken to be 0.05 « x\ « 0.1|| |j 
in the PQCD method. This suppression of the space-like penguins reflects the mechanism 
of the helicity-suppression, as stated in ||. 



To get the quantitative estimates, we take the wave functions as || 

M x ) = ^7^( Xl ~ ( 20 ) 

Mv) = V3/yJ/i(l-2/i), (21) 

M*) = Vzfzztil-zx). (22) 

Performing the integrations over [dx], [dy] and [dz], we get the final amplitudes: 



A ume = A(--21n- — ^ + ln^ — — - ztt) — 
2 2es £b £b 

3 1 

A time I A time A 

n 3d 3e — — o ' 

Z 6b 



A space = _ 936A) 

A space = 3A(1 _ ln l^_ m) 



(23) 



A time = 3 A(i-21n- — ^ + ln^ — f^-^) — , (24) 
2 2ts £b £b 



(25) 





BR(B~ -> n-K°) 


(xl(T 7 ) 


BR(B~ -> i^-iT ) 


(xlO" 8 ) 




■\xn "t" n /"in t" 

W 1 L11LJ Li L 


Txn"f" n 

W 1 L 11 


W 1 Lj 11LJ LI L 


ixn t" n 

W 1 Lll 




space-like 


space-like 


space-like 


space- like 




penguin 


penguin 


penguin 


penguin 


0.05 


8.38 


8.00 


10.3 


9.85 


0.06 


5.67 


5.36 


6.99 


6.61 


0.07 


4.10 


3.85 


5.05 


4.74 


0.08 


3.11 


2.90 


3.83 


3.57 


0.09 


2.45 


2.27 


3.02 


2.79 


0.10 


1.99 


1.83 


2.45 


2.25 



Table 1: Branching ratios for B — > tx K and B K K with space-like penguins 
included or not. 

when the diagram contributes to the specified decay process. The contributions of all the 

4^2 

other diagrams vanish. Here A = glfBfyfz- 

9 

The branching ratios for B~ — > K~K° and B~ — > vr^.ft' , are displayed in table 1. In the 
numerical calculations, we have used t b = 1.62 x 10~ 12 s, = 0.132GeV, = 0.113GeV, 
/ w = 0.093GeV, \V*V tb \ = 0.044 and \V t * d V tb \ = 0.013. 

From Table 1, it can be found that the space-like penguins contribute only several 
percents to the branching ratios and are thus less important. This result is consistent with 
the analysis we have carried out in the previous section. 

6 

We have carried out the analysis on the effects of the space-like penguin diagrams within the 
PQCD method. General suppressions of the space-like penguins are found compared to the 
time-like penguins. Numerically, the space-like penguins contribute only several percents 
in the penguin dominated processes B~ — > n~K° and B~ — > K~K°. We conclude that 



8 



space-like penguins are less important in this method. 
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Figure captions 

1. Time-like (a) and space-like (b) penguin diagrams. 

2. Graphs for the time-like penguins in PQCD. Here the dashed-lines denote the gluons, 
and the crosses denote the possible quark-gluon vertices. 

3. Graphs for the space-like penguins in PQCD. 
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